Journal of Nonlinear Mathematical Physics Volume 9, Number 1 (2002), |§-|75] 



Article 



Soliton Asymptotics of Rear Part of 
Non-Localized Solutions of the 
Kadomtsev-Petviashvili Equation 

Anne BOUTET de MONVEL t, Eugene KHRUSLOV^* , and Vladimir KOTLYAROV x * 

' Universite Paris-7, Physique mathematique et Geometrie, UFR de Mathematiques, 
case 7012, 2 place Jussieu, 75251 Paris Cedex 05, France 
E-mail: aboutet@math.jussieu.fr 

J Mathematical Division, Institute for Low Temperature Physics, 4 7 Lenin Avenue, 
310164 Kharkov, Ukraine 

* E-mail: khruslov@ilt.kharkov.ua 

* E-mail: vkotlyarov@ilt.kharkov.ua 

Received April 02, 2001; Revised November 03, 2001; Accepted November 04, 2001 

Abstract 

We construct non- localized, real global solutions of the Kadomtsev-Petviashvili-I equa- 
tion which vanish for x — ► — oo and study their large time asymptotic behavior. We 
prove that such solutions eject (for t — > oo) a train of curved asymptotic solitons which 
move behind the basic wave packet. 

1 Introduction 

The study of long-time asymptotic behavior of solutions of nonlinear evolution equations 
has attracted growing attention in last years and a number of papers have been devoted 
to this problem (see review paper || and references therein). Interest in this problem 
was especially stimulated by the discovery of the inverse scattering transform method 
[11, [l^, [l], |10| . In particular, one remarkable result obtained by this method was the proof 



that any localized (i.e. rapidly decreasing as |x| — > oo) solution u(x,t) of the Korteweg-de 
Vries (KdV) equation 

ut — 6uu x + u xxx = 



splits into a finite number of solitons when time tends to infinity [17]] (see also [18, jig]). 
Other nonlinear evolution equations with one space variable integrable by the inverse 
scattering transform also exhibit similar splitting. This phenomenon is an argument in 
favor of the physical interpretation of solitons as stable "long-living" particles. 

The simplest non-localized solution is of step-like form, i.e. the solution with following 
asymptotic behavior: 

0, x — > +oo 
— c 2 , x — * — oo. 



u(x, t) 
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It was proved in |l2| that the step-like solution of the KdV equation splits into [2^] 
soliton-like objects in the neighbourhood 

G+(t) = {x e M I x > 4c 2 t - Vint} 

of solution front (iV £ Z + and [ • ] denotes integer part). The form of these objects is 
similar to ordinary soliton but their velocities depend on i. In contrast with ordinary 
solitons they are not exact solutions of the KdV equation, however they satisfy it with 
increasing accuracy when t —* 00. For this reason such objects are called "asymptotic 
solitons". The number of these asymptotic solitons increases to infinity when t — > 00 
if the observation domain in the neighbourhood of the the solution front is extended 
correspondently. 

The same phenomenon of generation of asymptotic solitons trains on the solution front 
takes place (under certain conditions) for non-localized solutions of more general form 
as well as for other nonlinear evolution equations integrable by the inverse scattering 



transform [13]. Roughly speaking this phenomenon can be considered as a manifestation 
of the fact that any non-localized initial data consists of an infinite number of solitons 
which are gradually ejected at the front (first is the most rapid of them). An important 
condition for the ejection is the existence of wide "living area" for solitons where they can 
propagate without collisions. In the case where the non- localized initial data vanish for 
x — > 00 this area is a positive beam (at, 00) for some suitable a > 0. 

Another type of non- localized solutions of the KdV equation which vanish as x — > — 00 
was considered in |7|. It has been proved that under some conditions trains of asymptotic 
solitons are formed on a tail of the solution as t — > 00. These solitons move to the right 
following behind the basic wave packet. Physically it can be treated as an ejection of 
the slower solitons from a non-localized initial perturbation (initial data u(x, 0)). This 
phenomenon takes place under certain conditions on the spectrum of the Schrodinger 
operator whose potential is the initial data u(x, 0). This condition is the existence of a 
gap between the continuous spectrum of multiplicity one, which is provided by a non-trivial 
asymptotic behavior of u(x, 0) as x —* 00, and the continuous spectrum of multiplicity two 
(positive half-axis). 

In [||, ||, [6|, [TI| similar problems were studied for the Kadomtsev-Petviashvili equations 
(KP-I and KP-II): 

9 ( 3 1 \ 3 2 / N 

~dx r + 2 UUx + 4 Uxxx ) + 4 a Uyy = ' ^ ' 

where a = i for KP-I and a = 1 for KP-II. These are equations with two spatial variables 
x, y integrable by inverse scattering transform. See also (3|, Q for the Johnson equation 
and for the modified Kadomtsev-Petviashvili-I equation (mKP-I). 

Following the Zakharov-Shabat scheme of the dressing method we look for solutions in 
the form: 

u(x,y,t) =2^-K ± (x,x,y,t), (1.2) 
where the function K^(z,x,y,t) is a solution of the Marchenko integral equation: 

/■±oo 

K ± (z,x,y,t)+F{z,x,y,t)± K ± {s,x,y,t)F(z,s,y,t) ds = (1.3) 

J X 
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viewed as an equation in z > x (z < x) with parameters x, y, t. The kernel F(z, x, y, t) in 
(|1.3| ) satisfies the system of linear differential equations: 

+ F zzz = 0, 
= 0. 

Choosing solutions of this system in an appropriate way and solving the integral equation 
( |1,3| ) (with the sign "+" or "— ") one can construct solutions of the KP equations which 
vanish as x — > +oo or x — > — oo. For the KP-I equation global real solutions vanishing as 
x — ► +oo were constructed in || (see also Q). It was proved that in the neighbourhoods 
of their fronts: 

x>C(Y)t--?—logt N , 
a(Y) 

these solutions are asymptotically represented as follows: 

- / r 1 - x ~* 

u(x, y , t) = 2« 2 (Y) f cosh [a(Y) ( x -C(Y)t + — log t n ~ 1 / 2 ) - log <p n (Y) 

n=l 

+ 0{r l / 2+£ ) (i.5) 

where N is any natural number, Y = y/t, a(Y), C(Y) and (p n (Y) are some positive 
functions. This means that the trains of curved solitons are formed in the neighbourhood of 
the solution front because the ridges of these solitons at time t are located along curves x = 
C(y/t)t — a~ 1 (y/t) log i n_1 / 2 + log (p n (y/t) (n = 1,2,...). The phenomenon of formation 
of curved solitons trains nearly solution fronts is observed also for the KP-II equation, 
however in this case even solitons have a singularity ||. 

This paper is devoted to the study of non-localized solutions of the KP-I equation 
which vanish as x — > — oo. We construct a global real solution of this type and prove 
that a back part of the solution splits into curved asymptotic solitons. The dressing 
method of Zakharov-Shabat described above is not suitable for the investigation of long- 
time asymptotics of the tail of solutions. Therefore we use another method based on an 
integral equation in the plane of spectral parameters. This approach is also suitable for 
the study of long-time asymptotics of non-localized solutions in the neighbourhood of the 
solution front, which was first shown in |L4|] for the KdV equation. 

The paper is organized as follows. In Section 2 we prove the existence of a global real 
solution of the KP-I equation, vanishing as x — ► — oo. In Section 3 we reduce the problem 
to a degenerated integral equation and obtain a determinant formula for the solution. In 
Section 4 we study the asymptotic behavior of the determinant formula as t — > oo and 
prove an ejection of curved solitons which move behind the basic wave packet. 



2 Construction of global solutions vanishing as x — > — oo 

Let fi be a domain in the complex plane C (k = p + iq) with smooth boundary T = <9$7 
located in the right half plane at positive distance from the imaginary axis. Let us define 
the function E(p, q) := E(p, q, x, y, t) on $7 as follows: 



E(p,q) f{ P> q > Y W, 



(2.1) 
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with 



f(p,q,Y)=p 2 -3q 2 -2qY, 



(2.2) 



where x,y,t G R, and Y = y/t are parameters. Let us now consider the integral equation 
(with respect to ip(p, q) := ip(p, q, x, y, t)): 



4>(p,q) + 



1 



nX + k 



E(p, q)E{n, v)ip{n, u)g(fi, v) dfi dv = E(p, q), 



(2.3) 



where A = jU + wSfl, k = p + iq £ £1, k = p — iq and g(fi, v) is a smooth positive function 
on 0. Furthermore, if Q is unbounded, g(fi,v) satisfies 



^' + ^g(ji,u)dndu < oo 



(2.4) 



for any c > 0. 



Lemma 2.1. There exists a unique solution ip{PiQ.) = ^(PiQiZiUit) °.f (fe-^D which is C c 
with respect to x,y,t G R. 

Proof. denotes the Hilbert space of complex valued functions on Q with norm 



Let A be an operator on L 2 (Q), depending on parameters x,y,t G R, as follows: 



[A<p](p, q) 



E(p,q)E(p,v) 



n 



X + k 



^^)g(ix,v)d^ dv. 



(2.5) 



According to fl2.1|) , ( p.2j ) and (2.4), A is Hilbert-Schmidt for any values of the parameters 

-2t 



x,y,t. For ip G L g {Vt) we can write 



(A<p,tp) = [ E(p,q) [ 
Jn Jn 



X + k 



fiP, v)g(n, v) d\x dv <p(p, q)g(p, q) dp dq, 



where ( • , • ) denotes the inner product in L 2 {Q) and A = \i + iv, k = p — iq. Using the 
equality 



o 



X + k 

which is true because Re A > 0, Re k > 0, we obtain 



( x+ ' k >ds, 



(A<p,<p) 



ds 



n 



E(fi, v)(p(fi, v)e {ll+iv)s g(n, v) dfi dv 



> 0. 



(2.6) 



It follows from ( |2.6| ) that the operator A is positive. Therefore, the homogeneous equation 
ip+Aip = has only the trivial solution tp = 0. Since A is Hilbert-Schmidt, hence compact 
in L 2 (Q), the inhomogeneous equation 



V> + Atp = C 



(2.7) 
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has a unique solution in L 2 g (£l) for any ( £ L 2 g (Q). Let us note now that the integral 



equation (|2.3|) has the same form as ( |2.7[ ) with ( = E(p, q) = E(p, q, x, y, t) belonging to 
the space L 2 g (n) due to pip . Therefore, Q has a unique solution ip(p, q), depending on 
the parameters x, y, t. The first derivatives Dip of this solution with respect to x, y, t also 
satisfy (|2.7|) with right-hand side 



( = DE(p,q)- ! D 
Jn 



E(p,q)E{(i,u) 
X + k 



belonging to the space L 2 (^l). This proves their existence. Existence of high order deriva- 
tives is proved by induction. ■ 

Corollary 2.1. The inverse operator (1 + A) -1 exists and its norm in C(L 2 (Q)) is uni- 
formly bounded with respect to x,y,t G M: 



(i + A)- 1 !! < l. 



(2i 



Let us now define the function 



d f 

u(x, y, t) = -2— E(p, q, x, y, t)ip{p, q, x, y, t)g(p, q) dp dq, (2.9) 

where \/j(p,q,x,y,t) is the solution of the integral equation (|2.3| ), and E(p,q,x,y,t) is 
determined by (2J) and (^). 



Lemma 2.2. The function ( plj 
solution, vanishing as x — > — oo. 



is a solution of the KP-I equation (LI). It is a real 



Proof. Let us multiply (|2.3|) by E(p, q, z, y, t)e tq<yZ x ^g(p, q) and integrate with respect to 
p, q £ $7. Then, using 



1 



X + k 

together with (|2.1| ), (|2.2| ) we obtain 



where 



e (X+k)(s-x) 



i^(z, x, y,t) + / F(z, s, y, t)K(s, x, y, t)ds + F(z, x, y, t) = 0, 



(2.10) 



K(z,x,y,t) 



E(p, q, z, y, t)e %qi - z x) ip{p, q, x, y, t)g(p, q) dp dq, 



F(z,x,y,t)= / E(p,q,z,y,t)E(p,q,x,y,t)e iq{z x) g(p,q) dpdq. 
Jn 

Taking into accoun t (|2.1| ), ( |2.2| ), and (2.4) it is easy to show that the function F(z,x,y,t) 
satisfies equations (|1.4| ). It is also obvious that formula (2^) and equation fl2.10| ) (with 
respect to K(z, x, y, t), z < x) correspond to the equations (L2), Ql.3| ) (with the "— " sign). 
Thus according to the Zakharov-Shabat dressing method, u defined by ( [2,9| ) is a solution 
of the KP-I equation. Taking into account that the domain $7 is contained in the right half 
plane at positive distance from the imaginary axis and that the function g(p, q) satisfies 



inequality (2.4) it is easy to prove that solution (|2.9|) vanishes for x — > — oo. 
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This solution is real. Indeed according to (|2.9|) it is sufficient to prove Im(E, ip) = 0, 
where ifi is solution of the integral equation (|2.3|) and ( • , • ) denotes the inner product in 



Lg(£Y). Let us remind that equation ( |2.3| ) in L 2 (Q) takes the form ( |2.7| ) with A a positive 
operator and real right hand side £ = E. Then we obtain 



(E, = (C, ^) = (if> + A^, VO = W>, + (AV, ^) > 0. 
Hence Im(E, ijj) = and the lemma is proved. 



Thus for a given function g(p,q) > 0, formulas ( |2.9[) and (2,1)- (|2.3|) represent a global 
real solution u(x,y,t) of the KP-I equation. This solution decays as x — > — oo, but its 
behavior as cc — > +oo is unknown. 

Our main goal is to study the asymptotic behavior of this solution as t — » oo in suitable 
neighbourhoods of the rear part of the solution, namely neighbourhoods of the form: 

D N (t) = [(x,y) | -oo<y<oo, x < C(Y)t + log**}. 

We show that in such domains and for large t (t > T(N)) the solution u(x,y,t) has the 
asymptotic behavior described by (|L5|), where a(y), C(Y), tp n (Y) are expressed in terms 
of g(p,q). This formula means that the solution splits into a sequence of curved solitons, 
which are formed in the neighbourhood of its trailing edge. To prove this asymptotic 
formula we first need to approximate the solution of the integral equation ( |2.3| ) by solutions 
of integral equations with appropriate degenerate kernels. 

3 Integral equation with degenerate kernel 

Let ko = po + iqo G T be an arbitrary point of the boundary of $7. We will use the following 
double power series expansion: 

-. oo 

^=J2c ij (\-k ) i (k-k y, (3.i) 

i,j=0 

where 



i\j\{2 Po y+i+ 1 

It is easy to check that ( |3.1| ) converges in the polydisk 

n ( fc o) = {(A, k) | |A - feo| < po, \k- k \ < p }. 

Below we will choose ko as a point of f2 where the function f(p, q, Y) = p 2 — 3q 2 — 2qY 
attains its minimal value. We suppose such a point exists and is unique. This point 
depends on parameter Y: 

k = k (Y)=p (Y)+iq (Y). (3.2) 
Let us denote by C(Y) the value of the function f(p,q,Y) at the point ko(Y), i.e. 

C(Y) = f(p (Y),q (Y),Y)= min f(p,q,Y) (3.3) 
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and by Xn,y(j>, q) the characteristic function of a subdomain Gn,y C such that 

< dist(fco(n U \ G N , Y ) < ^p-. (3.4) 

This subdomain depends on Y and on N. It will be precisely defined later. 
Using the expansion (|3.l| ) we represent the operator A ( |2.5| ) in the form: 



A = An + + C]y + C 2 N , (3.5) 
where the operators A at, Bjv, C]y, and C 2 N are defined by 

[A N cp](p,q)= / E(p,q)E(iJ,,v)xN,Y(p,q)XN,Y{lJ>,v)g(fJ',v) 



n 

N 

x C{j(\ — ko) l (k — koyip(fi,v) dfidv, 
i,j=0 



^NV > ](p,q)= / E(p,q)E(n,v)xN,Y{p,q)XN,Y{n,v)g(ii,v) 



n 

x Cij(A — koY(k — /co)- J V(/ u i z/ ) dfidv, 

[CN<p){p,q)= / — — =-!— (1 - XN,Y(l*,v))g{ii,v)ip(ii,v)diidv, 
Jn \ + k 

[C%<p](p,q)= I E^l^^L^i _ XNY {p,q)) XNY (ii,v)g(n,v)ip(ii,v)diJ,dv. 
Jn \ + k 

Here A = fj, + iv, k = p + iq and 

RW := | i,j > 0} \ | < i,j < N} 

= | * > 0, j > N+ 1} U {(*, i) | t > N+ 1, i > 0}. 

Let us estimate the norm of the operators Bjv, Cjy and in the space Lg(Q). To avoid 
unessential complications we assume that Q is bounded. We will denote 

a = inf Re A, b = sup Re A, 

xen 

d(0 = (b + a)£ + (b - a)|£|, e = x - C(K)t, 

where C(Y) is defined by ( |3.3|) . The norms of the operators C l N (i = 1,2) on i^(O) can 
be estimated as follows: 

||C5v|| 2 < ff {\\ + k\)~ 2 e 2(p+ M e MC{Y)-f{p,q,Y))t e 2v,{C{Y)-f(w,Y))t 



x (1 ~XN,y(^ u )) 9{p,q)g{^,v)dpdqdndv 

e 2d(£)-2am t 



where 



^ (2a)2 $ (meas0) ' 



5= max ^(p, g) and mo = min [/(p, g, Y) — C(Y)]. 

(p,g)en ( P ,(j)en\Gjv,y 
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According to (|2.2f) , (|3J), (|3j), m > 0. Hence, 

IIC^II < 5 x measfi x e- amo * /2 (3.6) 

for £ < 2^*. Taking into account that C U fc >jv + i{(i, j) \i + j = k, i,j > 0} we can 
write (for A, k £ Gn,y)'- 

\Cij\ ■ \\ - koflk -kol 1 < ICijl-lX-ko^k-ko^ 

v ' i+j>N+l 

co I ^| 

EV — ttttIA - fcnl* \k - k \ l ~ l 
^ -i 2p ' +1 



1 oo 

— V 



I A — /cq I + I — I x 



< 



^ u l=N+l v ^ u 
1 /|A-/t | + |A:-A:o| x N+1 



Po V 2 Po 



Using this inequality we obtain 



where 



|Bjv|| 2 < (j e* \ ( — — —J xdpdqdfidv 

JJnxn Po V 2 Po ) 

= \ E iCijlJi&Y^Jjix^t), (3.7) 

P ° i,j>0, 
i+j=2(N+l) 



$ = 2(p + - t(/(p, g , y) - /(/x, f , y))] 
x = xjv,k(p, q)XN,Y{p-, v)g(p, q)g(n, v) 

Ji(x,Y,t)= I e 2 ^ x -^ p ' q ' Y ^\k-k \ i g(p,q)dpdq, (3.8) 



and the numbers Cjj are those introduced in (3.1). 

We now define more precisely the subdomain Gn,y C We will suppose that V is 
defined by 

r = {{p,q) I <p(p,q) = o} 

where ip(k) = <p(p,q) G C 2 (f2) and that the curvature of T is everywhere positive. Since 
the hyperbola 

H(Y):={(p,q)\f(p,q,Y)-C(Y)=0} 

is tangent to T at the point ko = ko(Y) E <9f2 we can introduce in some neighbourhood of 
ko new coordinates 

r:=2p(f(p,q,Y)-C(Y)) =: F(p,q) 

s := (|^)(p-po) " ^( fe o)(?-go))||V^fco)||- 1 =: *(p,g). 
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It is evident that s is the projection of (jp — po, q — qo) on the tangent to the boundary T 
at ko . It is easy to check that the equation of T near ko £ T is of the form 

r = a s 2 + 0(s 3 ), 

where 

ao = ||VF(fc )||^p. (3.9) 

ko and Ko are the curvatures of dtl and H(Y) at the point ko respectively. The minus sign 
occurs when T and H(Y) are on the same side of their common tangent. The plus sign 
occurs otherwise. 

In any case ao > 0. Hence, T is given by 

s = s±(r), r > 0, 

s±(r) = ±v^7ao + °( r )- 

Now, the subdomain Gn,y C O is defined as follows (Fig. [I]): 



Q 

Qo(Y) 


H{Y) 

1 _. r = 

k ^Xl^z — 




= e (N) 

r 






xV\\xv^^t— r 

Gn,y V 






Po(Y) j 





Figure 1. f2, T and the subdomain Gn,y 



G n ,y = {{r,s) | < r < e (N), s_(r) < s < s+(r)} 
where eo(N) is a small positive number, such that 

2 

£ o(^) < 7775 "°2° ,/> r27- ( 3 ' 10 ) 

16iV 2 (||(^o)| 2 + ||(fco)| 2 ) 



For £o(A0 > small enough both bounds in ( |3.4| ) are clearly fulfilled. 

Let us pass to an estimation of the integrals J{(x, Y, t) (i = 1, 2). Since for k G Gj\r y 



fc " fcol = 1/S 2 + „ V7rVI , 7JT2 (1 + 0( r + M)) 



||VF(fc 



o 
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we have 



Ik-ko] 1 < 2 l \s\ l + 



||VF(fc )||' 



Due to (|3.q ) we obtain 

Ji{x,Y,t) < 2 l e 2po ^ [ 
Jo 



eo(JV) 



e~ rt dr 



s+{r) 



s-(r) 



s l + 



\\VF(k W 



g(r,s)\w(r,s)\ds, (3.11) 



I) Op. q*\ 

where £ = x — C(Y)t, g(r,s) = g(p(r,s),q(r,s)), and w(r, s) = — ' ; is the Wronskian. 



Using Taylor's series one can write: 



d(r, s) 



e 2 ^ = l + 2(|(M + |(M)^o(r, S ,0 
where Eo(r,s,£) is bounded: 

\E (r,s,O\<\C\e 2 ^ Nm 
for (r, s) £ Gn,y, small Eq(N) and 



S (N) := 2||Vp(fc ) 



' e Q (N) 
a 



(3.12) 



(3.13) 



(3.14) 



Using this equality and the relations 

g(k) = g(ko) + 0(r + |s|) w(k) = w(k ) +0(r + \s\) 
we obtain after integration of ( |3.11 ) 



over s: 



Ji(x,Y,t) < — 



4 \ 2 ffoFol 



ao / i + 1 



eo(JV) 



e- rt r^(l + 0(v^))dr 



+ iei^(y)e 2po§+25oWI51 



i+2 



e r 2 



(l + 0(v^))dr 



<K 1J (y)- Tr + ^ 2i (y) 



i+4 



(3.15) 



where Ku and i^2i do not depend on £, i, and go = g(ko), wo = w(ko). Here we have used 
the estimate 



e- rt r d dr< I e~ rt r a dr 



-rtd r 







r(d + i) 



(3.16) 



where T(d) is the Euler T-function. From ( p.Tj ) and ( 3.1 5| ) we derive the following bound 
for the norm of Bjy: 



|Bjv|| <K N (Y) 



e 2p « e 2po^+S (N)\^\f 1 + \^\\ e 2po€+25 (JV)|e]n _|_ 1^1)2 



JV+4 
t 2 



JV+4 i 1 
i 2 t~4 



iV+5 

i 2 
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Therefore, in view of (gig ) and ( |1J) , we find 



l B ^ll ^f/S (0<e<l/4). (3.17) 



Now let us return to equation (|2.3|) (or (|2.7|) ). Taking into account (|2.5|) and fl3.5| ) we 
rewrite it in the form: 

ip = A N iP + B' N 4> = e (3.18) 
where = Bjy + Cjy- + and e = <?). According to (|3.6|) and ( 3.17| ) the norm of 



B^f has the estimate 



l^l^fi^r if t = x-C(Y)t<±-\ogt^. (3.19) 



Let us look for the solution of equation ( 3.18| ) in the form 

i> = Vjv + Sn, (3.20) 
where f/V is the solution of the equation 

iPn + AjvVjv = e, (3.21) 

and therefore 

5 JV = -(I + A)- 1 B^VJV- (3.22) 
According to Ql.2| ) and ( |3.20| ) the solution of the KP-I equation is represented in the form: 

u(x, t) = -2-^-(V, e) = -2-^-(Viv, e) - 2-^-(<5jv, e). (3.23) 
ax ax ax 

Taking into account (|3.22j ), the fact that A is self-adjoint, and relations (p. 18 ), ( 3.21[) we 

can write 

(5 N , e) = ((I + A^B'^iV, e) = (B'^iv, (I + A^e) 
= (B' N ip N ,ip) = (B' N ip N ,ip N ) + (B' N 7p N ,5 N ) 
= (B' N ^ N , i/> N ) - {B' N i; N , (I + A^B'^). (3.24) 

It follows from ( ^2l|) that 

IIV'Afll 2 + (ipN, A^at) - (lpN,B' N 1p N ) = (ipN,e). 
Hence, due to the positivity of the operator A, 

II^tvII 2 - \\B' N \\ ■ HVtvII 2 < \&N,e)\. (3.25) 
In the next section we will show that 

\& N ,e)\<KN, (3.26) 

where the constant K does not depend on x, y, t if £ = x — C(Y)t < 2^1ogt^^ +£ . 
According to ( 3.1 9| ) we have ||B^y|| — > as t — ► oo. Therefore taking into account ( p. 25 ) 
and ( pT26D we get 

\\ip N \\ 2 <2KN (3.27) 
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for sufficiently large t. In turn, it follows from ( |3.24D , ( |3.27| ) and ( j2.8|) that 

\(S N ,e)\<^l, (3.28) 

where does not depend on x, y, t if £ = x — C(Y)t < i^\ogt~^~ +e . 

Remark 3.1. The same estimate is also valid for the derivative ^(^AT,e). It can be 
proved using an analytic continuation of (ipN,e) in some strip | Imx| < [3. 

Thus according to ( |3.23D we need to investigate the solution of ( 3.21| ), This equation 
is an integral equation with degenerate kernel: 

N 

A N (p,q,fi,v,x,Y,t) = Cij(X - k y(k - k ) j 

i,j=0 

x E(p, q)E(n, v)xn,y(j>, q)xn,y{^, v)g{n, v) (3.29) 

and right-hand side 

E(p, q) = E{p, q, x, Y, t) = e P(*-/(P.* y )*) . (3.30) 
Due to the specific form of the kernel we look for a solution of ( |3.21 ) in the form: 

N 

i/> N (p,q,x,Y,t) =5^f°(s,y;t)(fe -ko) j E(p,q)xN,Y(p,q). (3.31) 

3=0 



Substituting ( |3.3l| ) into (|3.2l| ) and taking into account (|3.29|) , fl3.30|) we obtain a system 

•p-<r 



of linear algebraic equations for the functions = ip\ ■ (x,Y,t): 



N 

4 N) +J2 A jf^i N) =^ J=0,1,...,N, (3.32) 



1=0 

where Joo = 1 and Sjo = for j = 1, 2, . . . , N, 

N 



4 = 4jfa Y >t) =J2 c « J ^y^ ( 3 - 33 ) 



1=0 

and the integrals Jij are defined by 



J lj (x,Y,t)= E 2 (p,q)(k-koY(k-k o yg(p,q)dpdq. (3.34) 



The solution of the system ( 3.32j ) is given by 



ran d\ n \x, Y, t) 
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where (x, Y, t) = detf/^ +AS^ (x, Y, t)] is the determinant of the matrix with entries 
5ij + Aj!p(x, Y, t) (i,j = 0, 1, . . . , N), and D^ N \x, Y, t) is the determinant of the matrix 



obtained by replacing the j-th column of l( N ) + A^ by the column (1,0,..., 0) 



It follows from ( |3.31| ) and ( ]3.35|) that 



where F^ N \x,Y, t) is the determinant of the matrix obtained by replacing the first line 
of +A( N \x,Y,t) by the line (J (x, Y, t), J t (x, Y, t), . . . , J N (x, Y, t)). Let us now note 
that according to (|3.33| ), A^ N \x,Y,t) is the product of two (N + 1) x (N + 1) matrices: 

A^ N \x,Y,t) = C^j( N \x,Y,t). 

C (Ar) and j( N \x, Y, t) are the (N + 1) x (N + 1) matrices with entries Cy and Jy(x, Y, t), 
respectively. Taking this into account and setting Coo as a varying parameter in the matrix 
we obtain the determinant formula 

&N,e) = ^-logdet[jW +AW(x,Y,t)}. (3.37) 

4 Asymptotic behavior of the solution for large time 

First of all let us study the asymptotic behavior of the integrals ( |3.34| ). 

Lemma 4.1. The integrals Jij(x,Y,t) have the following asymptotic representation: 

an\wn\ h oH r ( ll± 2 ± ) e 2 ^ I ■(£ Y t) 

Mx,Y,t) = g ° KI ° i J [1 + (-1)^] 4rjTT + • 



Here, 



go = g(k ), w = w(k ), 

. I dk I f dp . dq\ 

ho = Tr( fc 0j = I Tj I" 1 ~7T~ ) (fco)> 

ao as ao \os as J 

where ao is defined in (|3ii|), and the functions Iij(£,Y,t) satisfy 

\Iij&Y,t)\ < Kij(Y)(l + |^|) e 2pof+MoWlCI 
with 5 (N) defined by ( gig ). 

Proof. Using ( 3,12| ) and taking into account that for A; € G/v,y 

k-k = ^(ko)s + |^( fc o)^ + 0(r 2 + |s| 2 ) 
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we write down the integral ( 3.34 ) in the form: 

f-eo(N) rs+(r) 



J ij (x,Y,t)=go\wo\k\k{e 2p °t e~ rt s i+j [l + 0{r\sf)]ds dr 

JO Js-(r) 

+ 2g \w \k[k{e 2p ^ e~ rt \p lS i+j+l + p 2 rs l+ ^ [1 + 0(r\sf)} ds dr 

JO Js^(r) 



s-(r) 

jh(x,Y,t) + 4(x,Y,t), (4.1) 



dk dp c)p 

where k\ = -z-(k ), Pi = Tr(ko), Pi = tt(^o) and [3 = if i + j = and (3 = — 1 if 
os as or 



i + j > 1 . Integration over s of the first summand in ( |4.1| ) gives the following asymptotic 
equality: 

T i ( v+ s go\w \k\k{e 2p °t i±i±i ^i+j+iin , ^ ni j 

J i -(x,y,t) = - ■ — 1 / e r 2 1 + 1 + O(Vr) dr. 

(i + j + l)a +J+L Jo 
Using the asymptotic relation ( |3.16| ) we find 



4(x,y,t)= ,. l w ! l 7, ~ - [i + (-ir w ]r 



(i + j + l)a \a J \a J L ' ' v 2 ; t ^±jt2 

In the same way, taking into account inequality ( |3.13| ), we obtain the following estimate 
for the second summand in ( |4.1| ): 

I4(x,y,t)| < KUY)\Z\ e 2po5+25o(7V)l5l / e -^ r a [1 + 0(^)}dr 

Jo 

<«y(y)^ 7+7+4 , (4-3) 

£ 2 

where the functions Kij(Y) do not depend on £ and i. The statement of the lemma follows 
from flOD-(p|). ■ 



Now, using ( |3.37| ), let us study the large time asymptotic behavior of the function 
#) (x, Y, t) = det[/M + ( Sj y, t)}. 



Lemma 4.2. VFe have the following asymptotic relation 



N 



2poC 



det[I (N) + A (iV) (X; y; t) ] = i + £ det det r<») (F)-^ [1 + * n (£, Y, t)} 

n=l t 2 

where and r^ n '(y) are n x n matrices with entries (i,j = 0, 1, . . . , n — 1) 



(i + jV- , 9o 

and 



i\j\{2 VQ yW ao (Y) i + j + l 

respectively, and the functions 5 n (£,Y,t) satisfy 

Kn(Y) ... 1 



6 n (Z,Y,t)\ ? /£<— logt^ +e with0<e<l/4. (4.4) 

tV4 2p 
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Proof. Let us denote by D^ N \x, Y, t; Ao, • • • , A/v) the determinant of the matrix + 
A^ N \x, Y, t), where A^) = diag(Ao, • • • , A/v) is the diagonal matrix depending on JV + 1 
parameters Ao, • • • , A/v- Clearly, 

D {N) (x, Y, t; 1, . . . , 1) = D (JV) (x, Y, t) = det \l {N) + ^ (x, Y, t) . 

This determinant is a polynomial with respect to the Afe's: 

W (x, Y, t; A , . . . , Ajv) = A . . . X N + A Ai . . . X N D^ (x, Y, t) 

+ A AiA 2 • • • X N D? (x, Y, t) + • • • + A . . . Ajv-i X N D ( ^ (x, Y, t) 

+ A AiA 2 • • • \ N D$ (x,Y,t) + --- + X ...X h ...X in ... X N D^ ln (x, Y,t) + . . . 

'0...N 



+ D { Q N) N (x,Y,t), (4.5) 



where D^ 1 ' i (x,Y,t) is the determinant of the n x n matrix with entries Ai r i p {x, Y, t) 
(r,p = 1, . . . , n); the hat means that the corresponding summand is absent. Taking into 
account ( |3.33[ ) and using Lemma 4.1 we obtain 



D hlin ( X > ^ *) = det ^ det J(n) ( X ' ^ *) + d[n) ( X > ^ *) 



= detcWdetrWfr)-^ + 4 n) (x,y,t), 

where j( n \x, Y, t) is the nxn matrix with entries J^\x, Y, t) defined by ( p.34| ), C^ n ^ and 
r( n ^(y) are defined in Lemma 4.2 (detC^ > 0, detr( n )(Y) > as Gram determinants), 
and the functions S n \x,Y,t), (x,Y,t) satisfy 



\dW(x,Y,t)\, \4 n \x,Y,t)\ < K d {Y) ^ 2) - [ 1 + ) . (4.6) 



The determinants 4 with ii + «2 + • • • + in > n ^" 2 ~ a ^ so satisfy (p~6|). Taking all this 
into account and setting Aj = 1 (i = 0, . . . , N) in ( [4.5|) , we obtain the assertion of Lemma 

L2|. ■ 

Remark 4.1. A more precise analysis of the determinants D^' in shows that the deriva- 
tives of the functions with respect to £ and to Coo have the same estimates as in ( |4.4[ ). 

Let us use the following equality, which is proved in Q: 

ndetcf = detc[ n ~ 1] , 

where is the nxn matrix with entries - ' ~rrir- (i, j = 0, . . . , n — 1), and c[ n ^ is the 

i. j. 

(i + J') ! 

(n — 1) x (n — 1) matrix with entries — — = 1, . . . , n — 1). This allows us to obtain 
the relation: 

^-detCH=2pondetCH| c ^ (2po) _ 1 , (4.7) 
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where is as above. 

Now taking into account ( 3.23| ), ( 13.28 ), ( 3.37| ), Lemma 4.2, Remarks 3.1 and 4.1 and 
I]) we obtain the following asymptotic formula for the solution: 



d 2 



i + J^detc (n) detr( n )(y)^ 



n=l 



(n + 2) 
t 2 



+ 0(V 1 / 4 ) (4. 



=x-C(Y)t 



in the domain {(x,y) £ R 2 \ x < C(Y)t + Jj- logi 2~ 3+e } as i — ► oo. 



2po 



Remark 4.2. This asymptotic formula is uniform with respect to y because, in ( 3.28p and 
(|4.4j) , Kn(Y) is uniformly bounded with respect to y. That follows from the compactness 
of the contour T and the positivity of its curvature. 



To push further the asymptotic analysis of the determinant formula (4.8) let us intro- 
duce notations: 



N 



A N (H,Y,t) = l + Y J Rn(Y) 



n=l 



R n {Y) = detC (n) detrW(y) 



£n(n+2)/2 



nil 



n-1 



(2 P0 ) 2n+2 n (*o 2 

i=0 



A ( n ) A ( n ) 
^1 ^2 i 



(4.9) 
(4.10) 



where Aj n) > 0, > are the determinants of the n x n matrices and with 

entries T(i + j + 1) and T((i + j + 3)/2) J> = 0, 1, . . . , n — 1), respectively. 

They are positive as Gram determinants. 
From (gl| and (O) it follows that 



a 2 a'^Aat - (a 

u{x,y,t) ~ ujv(a;,Y,t) = 2-—^ log A N (^,Y,t)\^ =x ^ C (Y)t = — 



I \2 

n) 



A 2 



(4.11) 



and 



AjvAtv - {A N ) - Ap ^ n(n+2)+i(i+2) • (4.12) 



n,i=0 



Let us cover the domain £ < ^ log i 2" +e by the intervals 



hit) 



-00 < £ < -^-logt 2+£ 

— log i 2 ~ e < £ < — log t 3+£ 
2p 2p 
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Taking into account (j4.1l| ) and ( |4.12|) , we obtain 



A 2 



R n ^ 1 (Y)e 2( - n - 1 ^ Rn(Y)e 2np °t 



(n-l)(n + l) 
t 2 



+ 



n(n+2) 



(i + o^- 1 / 2 )) 



and 



A'wA 



2R n (Y)R n . l{ Y)e^-^ f 



(Ajy) ^ ^ n ( n+ 2)+(n-l)(n+l) 



as £ G I n (i) and i — > oo. Hence, by virtue of ( p~9| ) and ( 4.11| ) we obtain 



u(x,y,t) = ^ 



2p 2 (Y) 



n= l COSh 2 

+ o(V 1/4 



P0 (Y) (x - C(Y)t + log^+V2 + x o (y) 



(4.13) 



where x° n (Y) = ^ log and F = y/t. 

Remark 4.3. The estimate (|3.26|) follows from (fD| and (p^ ). 

Thus we have proved the following result: 

Theorem 4.1. Zei i/ie contour F be compact, of class C 2 on £1, and with everywhere 
positive curvature. Assume the function f(p,q,Y) attains its minimal value at a unique 
point ko(Y) = po(Y) + iqo(Y) S T, for any Y = y/t. 



Then the solution u(x,y,t) of the KP-I equation defined everywhere by (|2.3j ) and (2.9) 
in the domain 

D N = {(x,y) | -oo < y < oo, x < C{Y)t + log t [("+D/2]+i+ e } (0 < £ < 1/4) 



has the asymptotic behavior defined by (4.8) and (|4.13[) /or £ 



oo. 
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